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Aba tract 

Calculations o£ the fan tone acoustic power and 
modal structuro genorated by complex distortions in 
axial Inflow velocity ate presented. The model used 
treats the rotor as a rotating thrcO-diincnslonal 
cascade and calculates the acoustic field from the 
distortion-produced dipole distribution on the 
blades including nonconpact source effects. Radial 
and circumferential distortion shapes arc synthe- 
sized from Fouricr-Dcsscl components representing 
individual distortion modes. The rclotlon between 
Individual distortion modes and the generated acous- 
tic modes is examined for particular dlatortlon 
cases. Comparisons between theoretical and experi- 
mental results for distortions produced by wakes 
from upstream radial rods show thot the analysis is 
a good predictor of acoustic power dependence on 
disturbance strength. 


Introduction 

Many papers have shown that the fan noise Isvel 
differences between static testing and projections 
to forward flight are due to the inflow distortion 
end turbulence present during most static test con- 
ditions,^ However, the importance of the Inflow 
distortion contribution to fan noise will depend on 
the particular fan stage and inflow distortion 
structure. It is therefore of interest to theoreti- 
cally study the fan noise due to inflow distortion 
interactions with a three-dimensional annular blade 
row and to compare the retiults of the tlieory to ex- 
isting experimental data. This theoretical study 
provides the fan noise modal structure which Is re- 
quired for the prediction of far-fleld radiation 
patterns^>3 and the design of acoustic suppressors 
where Inflow distortion Is the known source mech- 
anism. 4 

A previous paper, ^ presented the three- 
dimensional theoretical analysis of pure tone fan 
noise generated by Inflow Ulstortion/rotor Inter- 
action. It extended the thrcc-dimenslonel unsteady 
lifting surface theory developed by Nnmbo,® accepted 
as input various shapes of i.iflow distortion, and 
predicted the forward and aft radioted pure tone 
power and modal power distribution. Special empha- 
sis was placed upon the clarification of the role 
of the three-dimensional and noncompact source ef- 
fects. The numerical calculations for the full 
three-dimensional model, a two-dimensional model 
and a quasi three-dimensional model were carried 
out and compared for the case of an inflow distor- 
tion which was uniform radially and sinusoidal cir- 
cumferentially. However, real Inflow distortions 
are expected to have more complex shapes than the 
simple sinusoidal circumferential variation. There- 
fore, the present papaer explores the character of 
the predicted tone generation for more complex dla- 
tortlon types. In particular, the relation between 


the modal content of the Inflow distortion and the 
generated acoustic field is emphasized; and compari- 
sons between the theoretical and experimental re- 
sults for distortions produced by wakes from up- 
stream radial rods arc made. In the theoretical 
procedure, the local Inflow distortion is resolved 
into Fouricr-Dcsscl components, l.c. , inlet distor- 
tion modes. The effects of each inflow distortion 
mode are superimposed to obtain the individual acouE’~ 
tie modal powers generated which, In turn, arc 
summed to obtain the total pure tone acoustic power. 
An Important property of this Inflow distortion/ 
rotor interaction analysis is that the same acoustic 
mode can be generated by all the Inflow distortion 
components which have the same circumferential mode 
number, even though their rediel mode numbers are 
different. The next section Bumiuarlzes che theoret- 
Icel model which is described in more detail In 
Ref. 5. 


Analytical Model 

Fluctuating Velocity Induced on a Rotor Blade Row 

The theoretical model conslacs of a single 
three-dlmenslonol annular blade row with Njj blades 
rotating at constant angular velocity iii$ In an 
annular rigid-walled duct of infinite axial extent 
as shown in Fig. l(a), {Node that an .asterisk means 
the quantity is dimensional and lengths are non- 
dlmenslonslized by r$, the duct radius.) Thus the 
duct end reflection, the effect of upstream or do.wn 
stream blade rows, and the effects of duct area var- 
iation are not considered. But, these effects, ex~ 
cept for the area variation could be £' yzed using 
the solution procedure of the present utudy. The 
fluid flow (see Fig. 1(b)) is composed of an undis- 
turbed flow with a uniform axial velocity Hg and 
small fluctuating flows w^ ^ due to the inflow 
distortion. The flow is Invlscld, of uniform en- 
tropy and has no thermal conductivity. The fluctua- 
tions induced on the rotor by the inflow distortion 
convected by the mean flow are assumed to be Isen- 
troplc ond small compared with the undisturbed flow. 
It is also assumed that the fluid velocity relative 
to the blade is subsonic along the whole span and 
that the blades have no steady load, l.e.. a cascade 
of flat plate airfoils is analyzed. 

The flucuatlng velocity q Cr,8,z) induced 

by the rotor blade row/distortion Intaractlon can be 
obtained by integrating the linearized Euler's equa- 
tion of motion and the upwash component on the blade 
surface can be expressed in che following form: 


*NRC-NASA Lewis Research Center Associate; Natlonax Aerospace Laboratory, Tokyo, Japan, Member AIAA, 
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where Gq k ■X' acoustic dipole distribution on 
the blade* surfaces end KT|^ Is an upwash kcrnal 
function. The upwash kernal function contains pa- 
rameters of Hg (blade number), h (hub- tip ratio), 
q and p (circumferential and radial mode numbers 
of Inflow distortion), (H'p (rotor spccd/axlal flow 
speed), 0 (Interblade phase angle) and L (number 
of terms In the finite scries approximation for 
•'o k)‘ 7;he detailed expresfllonQ for the kernal 
function ECTk are given In Ref. 5. 


Inflow Distortion 


In this paper, the inflow distortion Is assumed 
to have only an axial velocity component. (Other 
components con be considered using the same theoret- 
ical methods.) When n Fourler-Bessel analysis of 
arbitrary shapes of inflow distortion Is carried 
out, the following axial component of external 
fluctuating velocity Is obtained; 


w^,^(r.e,z,t) 
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where R_(kqjpr) Is a radial eigenfunction of q 
order, Cq denotes a small quantity which Is the 
ratio of the external fluctuating velocity to mean 
flow velocity and Bq „ are the Fourier coeffi- 
cients of the Inflow distortion. Then, the upwash 
component of the external fluctuating velocity on a 
blade surface con be expressed by : 
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r, is convenient to suppose that the distortion 
velocity can be expressed as the product: 


then 
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When the Inflow distortion can be represented by H 
Gaussian profiles (for N 1, 2, ...) In the cir- 
cumferential direction, then! 


Oj^(S) - exp 


and we find that: 
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0 otherwise, for q “ 0, +N, +2N, ... 
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Model of Rod Wakes 

For the specific case of an axial velocity dis- 
tortion produced by the wake of a cylindrical tube 
or rod of diameter d*i the magnitude of the wake 
defect, Ej,, and the Gaussian half-width of the wake 
profile, 02 , are given by the following equations 
based on Prandtls’ mixing length hypothesis (cf., 
Ref. 7, p. 691). 
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where 0j^(6) is given by eq. 8 and 
a _ 0.298 .* * 

“ 2xri r-' ^ 


( 10 ) 


( 11 ) 


The N rods are a distance x* upstream of the 
rotor blade tip leading edge and Cg is the drag 
coefficient of a rod. 


Determination of Acoustic Dipole Distribution 

The upwash component of the external flucuatlng 
velocity q^ y must be cancelled by the Induced 
upwash velocity q^(r,0,z)c^^ at the blade sur- 
faces. Then, an integral equation for the unknown 
acoustic dipole Gg^k(^>q>p) obtained In the 
form: 
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0 ( 12 ) 


The quantity Gg „( 5 ,q,p) is aalculated from 
cq. (12) by u collocation muthodi 


Pure Tone Acoustic Poucr 


Thu dlmenslonlcaa acoustic power with ro- 
apect to the Jth pure tone fan hamonlc (for ex- 
ample, J « ±1 corrospondc to the fundamental) la 
given by: 


and 


4 



Ell(n,l,q) 


(13) 


El: 


j |HPj,(n,l,q)r B^(nw^ + m) J 

(a^ + nuy lit) ^ 


(W 


El^(n,f,q) Is the modal component of the dlmensluu- 
loB8 acoustic power, and is nondlmenslonallzed with 
s/4 pjjTT^c^rj^, llPj(n,)l,q) dent tea the nondlmen- 
alonal p'-essure amplitude In the (n,S.) acoustic 
mode, which is given by: 
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X expj-l(a^ + nuj^)cjdc ( 15 ) 

a± = Is the axial wave number nondlmenslonallzed by 
tj. In eqs, (13) to (15) the + and - eubscrlpts 
refer to downstream and upstream propagation, re- 
spectively. For Ihe compact source analysis, 
exp[-l(a+ + nlil'r)cj » 1 In eq. (15). 


The numerical calculations are carried out in 
this paper for three different rotor configurations. 
Their overall geometric and aerodynamic parameters 
arc given In Toble I and the radial variations of 
otaggor angle, y, pitch-chord ratio, S, and ratio 
of blade relative to axial velocity, q/T7g, ore s)ujwn 
In Fig, 2. The numerical results Include four cases; 
(1) the acoustic power In the modal components gen- 
erated by particular Fourlcr-Deasel components of 
inflow distortion In the case of the fop denoted 
No. 1; (2) the acoustic power and modal ,':ontent 
variations with the half-width of a slngl.t Gaussian 
profile representing the local Inflow distortion 
profile In the circumferential direction (fan 
Ho. 2); (3) the comparison between thcorjtlcol and 
experimental tone power for the case of on upstream 
distortion produced by the wake of a cylindrical 
tube Immersed to varying depths from the duct wall 
(fan 2) ; and (4) the comparison between numeri- 
cal nud experimental tone powers and the analysis 
of modal content for the cases of 28 and 41 rod 
wakes interacting with a 28 bladcd rotor (fon 
No. 3). 


Acousclc Modal Power Generated by Portlcular 
Fourler-nossel Inflow plstortlon Modes 


Figs. 3 and 4 show the fundamentol pure tone 
modal powers as a function of the circumferential 
lobe number q of the Inflow distortion whose 
radial distribution Is assumed to be the 1st 
(Fig. 3) or 3rd (Fig, 4) radial eigenfunction of 
order q. Figs. 3(a) and 4(a) are for upstream 
propogatlon, while Figs. 3(b) end 4(b) are for down- 
stream propagation. The fan used in this calcula- 
tion 1ms the largest number of blades of the three 
cases as shown In Table I and, for similar tip speed 
conditions, gcncr.ites the largest number of propa- 
gating modes. The calculations were limited to 
q £ 50, The 1st radial inflow distortion mode has 
a maxiraiira amplitude at the rotor tip, while the 3rd 
radial inflow distortion mode hos pcoks near the 
middle of the blade span. With the radial distor- 
tion mode number p held constant, the radial dis- 
tribution of Inflow distortion velocities changes 
continuously with circumferential distortion mode 
number q with more skewing toward the wall as q 
Increases, The maximum amplitudes of all Inflow 
distortion modes used in these calculations are 
equal to 1. The numerical results show that while 
the 1st radial mode of Inflow distortion generates 
all radial acoustic modes, the 1st radial acoustic 
mode, t X 0, dominates for most values of q. The 
3rd radial mode of Inflow distortion also generates 
nil radial acoustic modes; but, over a wide range of 
q, m.iny of them carry mote power than p « A n 2, 

In the case of the 1st radial mode of Inflow 
distortion (Fig. 3) , the 1st radial acoustic mode 
power is larger than other radial acoustic mode 
powers by more than 5 dB In both upstream and down- 
stream casus, except for the points at q = 40 and 
the downstream values at q = 50. At q 40, the 
circumferential Inflow distortion mode number equals 
the rotor blade number and the circumferential 
acoustic mode nvunVer, n, Is zero (n = jNg - q). The 
(0,0) acoustic wod.< I.' a ulare wave. If one accepts 
the Intuitive nut lot thid the acoustic pressures in 
the various mo us wlM s :m in a manner such that 


the net rodlol, verlotlen will be Blnllar to the 
radial varlotlon of dlatortloni then It lo rcoeon- 
able to Dxpoct a combination of 2nd and 3rd radial 
acoustic modes will exist at q 40 in addition to 
the plane wove. For Che cases where n 0, the 
i - 0 modes opparontly provide o reasonable radial 
match to the (qi 0 ) distortion modes. 

In the COSO of the 3rd radial mode of inflow 
distortion (Fig. 4), the acoustic powers of the 
higher order radiol modes (4 > 2 ) contribute to the 
total fundamcntol tone power over the midravgc 
values of q which correspond Co pcok total funda- 
mental power generation. The 3rd radial acoustic 
mode (t •• 2 ) dominates only at lower values of q 
near 20. Note that the t >• 2 acoustic mode in 
Fig. 4(a) shows a distinct minimum at n ■■ 10. This 
is probably a point where the resultant acoustic 
dlpola distribution on tlie blades cannot couple 
with the ( 10 , 2 ) acoustic mode similar to the case 
discussed in connection with Fig. 13 of Ref. 9 
which applied to rotor-stator interoctlon. In gen- 
eral, for cose' other than the sljtplest (p •• 0 ) 
radiol varlotlou of Inflow dlatortloni we arc not 
yat able to rationalize the calculated distribution 
of radial acoustic mode powers produced by a single 
higher order radial distortion mode. 


Acoustic Power and Hodnl Content Generated bv o 
Gnusalan Circumferential Inflow Distortion Profile 

Fan No. 2 used in this calculation has a low 
tip speed and small blade number. Therefore, the 
totol number of acoustic modes generated at a par- 
ticular circumferential distortion mode number is 
less than for fan Ho. 1. Fig. 5 shows the acoustic 
power generated by fan No. 2 intcraotlng with n sin- 
gle Inflow distortion represented by a Gaussian pro- 
file in the circumferential direction and a radial 
step function velocity defect that extends from the 
blade tip down to 72 percent of the span. The power 
is plotted as a function of the half-width of the 
Gausslon profile, 8 2 * Results are shown for the 
total, 1st, 2nd, and 3rd harmonic cone powers pro- 
pagating in both upstream and downstream directions. 

Fig, 5 shows that the acoustic power in any 
harmonic peaks at particular values of the Gaussian 
half-width. This phenomenon, mentioned in Ref. 8 
and discussed in Ref. 1, Is related to an interplay 
between the number of propagating acousCi: modes and 
the Fouriet-Bessel inflow distortion harmonic am- 
plitudes which can couple to the acoustic modes. An 
example of the circumferential Fourier distortion 
harmonics corresponding to the conditions of Fig, 5 
is shown in Fig. 6 for several values of Gaussian 
profile helf-wldth. The range of n (or q) corres- 
ponding to propagating circumferential acoustic mode 
numbers Is also shown. With an increase in 82 , the 
amplitudes of the lower order circumferential dis- 
tortion coefficients, |aq| 2 , increase; but the large 
amplitudes occur over a smollct portion of the pro- 
pagating range. The result of this coupling con- 
straint is a peaking and decline beyond a certain 
value of profile width. 

For Values of 82 greater than about 0.02, 
the downstream propagating acoustic power in the 
fundamental tone ia higher than the upstream power. 
For the 2nd and 3rd harmonic and the total power, 
the downstream power is olweys greater than the 
upstrerm power. Of course, in a real fan, down- 


stream radiated power must poss through a stator 
blade row and is subject to modification. 

Figs. 7(a) and (b) compare compact source pre- 
dictions with the noncompaci calculations of Fig. 5 
for upstream and downstream propagation, respective- 
ly. These figures indicate that the ocoustic power 
differences between corresponding harmonics are con- 
siderable for the fundamental and 3rd harmonic pro- 
pagating upstrenm and for all harmonics propagating 
downatream, partlculorly in the range of 62 around 
the peak acoustic power generation 

An explanation for this behavior is as follows: 
Uic results in Fig. 7 correspond to summations over 
the inflow distortion circumferential mode number 
q. As q is varied, the noncompact ocoustic power 
in the fundamental tone la sometimes higher end 
sometimes lower than that based upon the compact 
prediction (see Fig, 7 in Ref, 5). At large volues 
of 62 , the number of circumferential distortion 
harmonica which )n couple to propagating acoustic 
modes decreases and the net acoustic power differ- 
ences between noncompact end compact results arc 
considerable and vorled in magnitude. At small 
values of 62 , the number of 'Coupling distortion 
harmonics increasca and therefore, the acoustic 
power differences among modes ate averaged out to 
small not amounts in the summation over the wide 
range of coupled q's for o portlcular tone har- 
monic. The compact source predictions In Fig. 7 
tend to underestimate the upstream acoustic power 
ond overestimate downatvenm acoustic power in the 
fundamental tone. The trends at the 2nd ond 3rd 
harmonics indicate that 0 compact source assumption 
underestimates harmonic power for both propagation 
directions. These results differ from those of 
Ref, 9 which considered only trends associated with 
single sinusoidal modes of Inflow distortion enter- 
ing B two-dimensional cascade. The summation over 
many distortion modes representing more complex 
shapes such as the Gaussian, makes the effects of 
neglecting noncompactnesa specific to Che particular 
case considered. 

Figs. 8 (a) and (b) show the variation of the 
acoustic mode structure with the half-width of the 
Gausslon profile, 82 . Note that the minimum power 
calculated was -70 dU end that lower values are 
plotted at -70. The step radial tip distortion as- 
sumed favors 1 st radial mode generation similar to 
the results discussed in connection with Fig. 3. 

For particular modes , clear generation minima are 
also evident, o.g. , for Bg " 0.0125, 4 » 0 , n " 2 
in Fig. 8 (a). The envelopes of Che peak powers of 
the individual modes follow the trends with q 
shown in Fig. 6 , 


Cociparlsons of Theory and Experiment 
A Single Tube Wake InteractinR with e Rotor 

Experiments were conducted in an anechoic wind 
tunnel^® in which the variation of upstream todlated 
fundamental cone power was measured as function of 
the length of a cylindrical probe tube inserted 
through inlet wall upstream of rotor No. 2. The 
wind tunnel was operated with 40 knots velocity and 
the fan fundamental tone due to rotor-statot Inter- 
action was cutoff such that inlet tone noise was 
dominated by the probe tube wake interacting with 
the rotor. Fig. 9 shows the measured and calculated 
tone power (referenced to full immersion) as o func- 
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ticn of prob<^ Imnorelon. The detailed Bhnpc of tlic 
ptobo tip Is lllufltracod In bo upper right corner 
of Flgi 9* The probe tube wake was modeled with 
Eqa. (9) CO <11) as having a Gauaslon profllo in 
the circumferential dlrccclon. A step profile cor- 
responding to Immersion length was used to represent 
Che radlol distortion vorlotlon. The theoretical 
results show good agreement with the measured In- 
creases in fundamental cone power per length incre- 
ment of probe Immersion. For full immersion) I. “ 

11 cm, Che calculated tone power is 130 dB compared 
with the measured value of 133 dD. The absolute 
amplitude calculated is subject to some uncertainty 
associated with the wake modeling; c.g. , Che effec- 
tive O 2 used CO rcprcscnC the complex probe tip. 

Of more significance is the agreement in the mea- 
sured power trend with Incroablng probe Imraorslon. 

Fig, 10 shows the variation of acoustic mode 
powers contributing to the fundamental tone as a 
funcClcn of probe tube immorslon. The differences 
in mode power between 6.7 cm and 11 cm immersion 
are seen to bo small for both the 1st and 2nd radial 
modes echoing the overall power variation shown in 
Fig, 9 ond emphasizing that the generated power is 
dominated by the fan tip region as previously shown 
in Ref. 5. 


Multiple Rod Wakes IntcrnetlnR with a Rotor 

A JT15D turbofan engine was operaced on an out- 
door teat stand with two separate sets of distortion 
rods (28 and Al) as a means of studying the trans- 
mission characteristics of inflow control de- 
vices. The inflow control greatly reduced fan 
sources associated with inflow disturbonccs and left 
the dominant rod wake - rotor interaction os the 
main source of inlet noise. The fan parameters and 
operating condition arc described in Table I as fan 
Mo, 3. The 28 rods were 0,635 cm in diameter and 
extended from the wall 57,8 percent of the blade 
span, and the 41 rods were 0.476 cm in diameter and 
extended from, the wall to 62.5 percent span. The 
centerline of both rod sets was located 10.3 cm up- 
stream of the rotor tip. Fig, 11 compares the cal- 
culated and measured fundamental tone power and 
shows the calculated distribution of modal powers 
for the two cases. The modal powers ore roughly 
equal for the 28 rod case with the exception of the 
low 3rd radiol which appears to be another example 
of weak coupllcg between blade dipoles and the (0,2) 
acoustic mode. The agreement between measured and 
calculated tone powers is rather good in each case, 
but such ogreements of Isolated powers are subject 
to some skoptlclsra. The quantity of more Interest 
is the comparlscn of fat-field directivity between 
measurements and predictions from the calculated 
modal content. However, several links in the ana- 
lytical chain relating source modes in an annulus 
to far-field directivity are missing. The modal 
scattering in the transition from an annular to a 
circular duct awaits detailed numerical solution os 
does the detailed radiation from an unflanged inlet 
duct to the far-field with a superimposed inlet po- 
tential flow field. Ref. 12 looks at this same rod 
disturbance data from the standpoint of Inferring 
source modal content from the far-field directivity. 
A reconciliation of those results with the present 
calculation requires additional analysis. 


Concluding Remarks 

The three-dimensional, noncompact source theory 
has been applied to calculate fan tone power and de- 
tailed acoustic mode structure for several assumed 
inflow distortions of increasing complexity and for 
experimental situations whore upstream radial tods 
generated the distortions. While generated radial 
acoustic mode structure can be rationalized for tip 
distortions (lot radlnl dlscortion modes appear to 
favor 1st radiol acoustic modes), the radial acous- 
tic mode content becomes much more complex for high- 
er order radial distortion modes. Circumferential 
distortion profiles such as the Causslan used herein 
simply generate circumferential acoustic mode con- 
cent corresponding to the circumferential Fourier 
harmonic content of the distortion. However, indi- 
vidual acoustic mode levels are subject Co a cou- 
pling constraint associated with the way the dipoles 
on Che blade surfaces couple to a particular acous- 
tic mode. Poor coupling results in very low values 
of generated mode power. The use of this minimum 
gcnecoclon property to prsctlcally reduce fan noise 
would require that the mode minimized be the sole 
mode responsible for the particular tone harmonic 
power generation at the given fan operating condi- 
tion. The presence of other modes not obeying the 
particular decoupling constraint would defeat the 
tone power minimization. Generolizacions regarding 
Che magnitude and sign of the differences between 
compact and the noncompoct source results ore not 
available for complex Inflow distortions since the 
net tone powers are the result of summations over 
many individual modes, each having varied degrees 
of noncompact dependence. The fact that the anal- 
ysis is able to predict the trend in Cone power ns n 
function of the immersion depth of an upstream rod 
disturbance Is encouraging. An important task re- 
maining is to link acoustic mode content at the fan 
face to far-field directivity patterns which can be 
compared with experimental results. 
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TABhE I. - FAN PARAMETERS 


Fan number 

1 

2 

3 • 

Rotor blade number, Ng 

40 

15 

28 

Ptmcnslonless axial tip 

,055 

.258 

.094 

chord length, 




Rotor tip specd/axial 

2.474 

1.052 

3.011 

flow speed, 




Rotor relative Mach 

.934 

.865 

,913 

number, Mx 




Axial flow Mach number, 
Ma 

Hub/tlp ratio, h 

s350 

.596 

.288 

.40 

i 

.40 
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Figure 3. - Fundamental pure tone power variation with 
radial inflow distortion given by the Fourier-Bessel 
1st radial mode of order q. fan no. 1, maximum 
amplitude • 1. 
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Figure 4. - Fundamental pure tone power variation wItO 
radial inflow distortion given by the Fourier-Bessel 
3rd radial mode of order q, fan no. 1. maximum 
amplitude • 1. 
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Figure 7. - Comparison of compact and noncompact source 
acoustic power prediction lor a single Gaussian circum- 
ferential distortion, fan no. ?, 28% tip radial distortion. 
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Rgure 8. - Variation of the fundamental tone modal 
structure with half width of Gaussian circum- 
ferential distcrtion profile. 28% tip radial distor- 
tion. fan no. 2. 




Figure 9. - Fundamental tone acoustic power generated by 
rotor-probe tube wake interaction. fann >. ?, upstream 
propagation. 
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Figure 10. - Variation of fundamental tone modal structure with 
probe length immersed in a fan inlet, fan no. 2, upstream 
propagation. 
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Figure 11. - Modal structure of the fundamental pure 
tone power generated by rod wakes interacting with 
a fan, fan no. 3, upstream propagation. 



